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1. Introduction . ; . : : : 
devices is one of the most important technical issues in many 


high-tech industries. The conventional method to increase the 

cooling rate is to use extended heat transfer surfaces. How- 

ever, this approach requires an undesirable increase in the size 

of the thermal management system. In addition, the inherently 

poor thermal properties of traditional heat transfer fluids such 

ae oo as water, ethylene glycol (EG) or engine oil (EO) greatly limit 
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With an ever-increasing thermal load due to smaller fea- 
tures of microelectronic devices and more power output, cooling 
for maintaining desirable performance and durability of such 
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Nomenclature 
Cop an empirical coefficient in Rohsenow correlation 
Cp specific heat (J/kg K) 


D tube diameter (m) 

f friction factor 

g acceleration of gravity (m/s2) 

h heat transfer coefficient (W/m? K) 
k thermal conductivity (W/m K) 

K a constant in Zuber correlation 

L length (m) 

m mass flow rate (kg/s) 


q supplied heat (W) 
q" heat flux (W/m?) 
T temperature (K or °C) 


Greek symbols 


(o particle volume fraction 
density (kg/m?) 

H dynamic viscosity (kg/m s) 

o surface tension (N/m) 

Subscripts 

b bulk 

f fluid 

g vapor phase 

i inner 

o outer 

nf nanofluid 

p particle 

s saturated liquid 

v laminar sublayer 

w wall surface (heater or tube) 


Nondimensional numbers 


Nu Nusselt number 
Pe Peclet number 

Pr Prandtl number 
Re Reynolds number 
Abbreviations 


BHF burnout heat flux 


BHTC boiling heat transfer coefficient 
CHF critical heat flux 

CNT carbon nanotube 

DNB departure from nucleate boiling 
Ni-Cr nickel-chromium (nichrome) 
ONB onset of nucleate boiling 


a vital role in the development of energy-efficient heat trans- 
fer equipment it is known that fluids possess order-of-magnitude 
smaller thermal conductivity than metallic or nonmetallic particu- 
lates. For example, thermal conductivities of water (0.607 W/m K) 
and engine oil (0.145 W/mK) are about 5000 times and 21,000 
times, respectively smaller than that of carbon nanotubes (e.g. 
3000 W/mK for multi-walled carbon nanotubes, MWCNTSs). Fig. 1 
compares thermal conductivities of various commonly used liq- 
uids and materials at room temperature. It can be seen from Fig. 1 
that most of the commonly used metallic and nonmetallic par- 
ticulates posses orders magnitude higher thermal conductivity 
than those of liquids. Fig. 1 also shows that metals have higher 
thermal conductivity compared to ceramics (or metallic oxides). 
Therefore, the thermal conductivities of fluids that contain sus- 
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Fig. 1. Thermal conductivities of commonly used liquids and materials at room 
temperature (some small thermal conductivity values are shown above the numbers 
of the liquid or material) [1-6]. 


pended metallic or nonmetallic particles or tubes are expected 
to be significantly higher than those of traditional heat transfer 
fluids. 

Although nanoparticles-suspensions were used in heat trans- 
fer studies as early as 1984 by Yang and Maa [7] and then 
in 1993 by Masuda et al. [3], it is only in 1995 that Choi [8] 
at Argonne National Laboratory of USA coined the concept of 
“nanofluid” to meet the cooling challenges facing many high- 
tech industries. Apart from Masuda et al. [3], there are other 
researchers who also used the term “nanofluid” or were involved 
with research on thermal conductivity of nanoparticle suspen- 
sions during 1993. Gass et al. [9] from Switzerland used the 
same term “nanofluid” to express minute volume of fluid (nano- 
liter) in microfluidics study in 1993. A German researcher, Grimm 
[10] also won a German patent on the enhanced thermal con- 
ductivity of nano- and micro-sized particles suspensions in the 
same year. Aluminum particles of 80 nm to 1 um were suspended 
into a fluid and about 100% increase in the thermal conductiv- 
ity of the fluid for loadings of 0.5-10vol.% was reported in his 
patent. 

The definition of this new class of heat transfer fluids, i.e. 
nanofluids is the suspensions of nanometer-sized solid parti- 
cles, rods or tubes in traditional heat transfer fluids. Nanofluids 
were found to exhibit significantly higher thermophysical prop- 
erties, particularly thermal conductivity and thermal diffusivity 
than those of base fluids [11-16]. Nanofluids have attracted great 
interest from the research community due to their enhanced ther- 
mal performance, potential benefits and applications in numerous 
important fields such as microelectronics, microfluidics, trans- 
portation, manufacturing, medical, and so on. Recent record shows 
(Fig. 2) that there is an exponential growth of annual research 
publications on nanofluids. According to ISI Web of Knowledge 
searched results (Fig. 2), total 153 nanofluids-related publications 
have already appeared in the first half of year 2010 (till July 
14). It is also believed that there are more than 300 research 
groups and companies worldwide who are involved with nanoflu- 
ids research [17]. Compared to studies on thermal conductivity, 
few works have been reported on boiling, droplet spreading, and 
convective heat transfer characteristics of nanofluids although 
these features are very important in order to exploit nanofluids 
as the next generation coolants. It is, therefore imperative and 
timely to provide a state-of-the-art review on the advances in boil- 
ing, droplet spreading, and convective heat transfer research of 
nanofluids. 
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Fig. 2. Annual publications on nanofluids since 2000 (publications including all 
types of journal and conference articles, patent, news, letter and other). 
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Fig. 3. Typical pool boiling curve. 


2. Boiling and droplet spreading of nanofluids 
2.1. Boiling heat transfer—basic 


Boiling is a phase change (from liquid to vapor) process in which 
vapor bubbles are formed either on a heated surface or in a super- 
heated liquid layer adjacent to the heated surface. Although boiling 
is a complex and elusive process, it is very efficient mode of heat 
transfer in various energy conversion and heat exchange systems as 
well as cooling of high energy density electronic components. There 
are two common types of boiling: pool boiling and flow or forced 
convective boiling. Pool boiling refers to boiling on a heated surface 
(heater) submerged in a pool of initially quiescent liquid while flow 
boiling is boiling in a flowing stream of fluid, where the heating sur- 
face may be the channel wall confining the flow. Heat flux in any 
boiling is one of the most important parameters in designing and 
operating the heat transfer equipments with high heat flux such as 
boiler, evaporator, electronic equipment and rocket engines. Since 
there are very few efforts are made on flow boiling of nanofluids, 
this mode of boiling heat transfer will not be elaborated further 
instead pool boiling heat transfer is briefly discussed here. 

In pool boiling problems, the fluid is initially quiescent near 
the heating surface, and subsequent fluid motion arises from free 
convection and the circulation induced by bubble growth and 
detachment. Two main parameters in boiling problems are the 
degree of wall superheating or excess temperature (the difference 
between the wall and the liquid saturation temperature at the local 
pressure, AT = Tw — Ts) and heat flux (q”). The classical pool boil- 
ing curve, which is as a plot of heat flux versus excess temperature 
(AT), is shown in Fig. 3. As the value of AT increases, the curve 


traverses four different regimes or modes of pool boiling (Fig. 3): 
(I) free convection, (II) nucleate boiling, (III) transition boiling, and 
(IV) film boiling. 

As illustrated in Fig. 3, up to point A is the free convection regime 
which occurs since there is insufficient vapor to cause active boil- 
ing. Small temperature differences exist in the liquid, and heat is 
removed by free convection to the free surface. At point A, isolated 
bubbles initially appear along the heating surface and this is termed 
as onset of nucleate boiling (ONB). Nucleate boiling occurs between 
points A and C. At this regime, isolated bubbles appear, and heat is 
transferred mainly from the surface to the liquid. As AT increases 
(B-C), more nucleation sites become active and bubbles coalesce, 
mix, and ascend as merged jets or columns of vapor. At point C, 
the maximum heat flux or critical heat flux (CHF) occurs. Between 
points C and D, transition boiling occurs. An unstable (partial) vapor 
film forms on the heating surface, and conditions oscillate between 
nucleate and film boiling. Intermittent vapor formation blocks the 
liquid (higher conductivity) from contacting the surface (lower- 
ing the surface heat flux). Film boiling occurs beyond point D. In 
addition to conduction and convection, heat transfer by radiation 
is important at these high wall superheat levels. A stable vapor film 
covers the surface in this regime. 

There are numerous factors such as heater or channel surface 
conditions (smooth or rough), heater size, shape, material, diameter 
and orientation, degree of surface wetting, subcooling, inclusion of 
surfactants, and properties of liquid affecting heat transfer perfor- 
mance and bubble dynamics of pool and flow boiling. These factors 
are widely studied both theoretically and experimentally [18-23]. 
For instance, effects of several parameters such as tube diame- 
ter, length, surface roughness, and inclination on pool boiling heat 
transfer have been investigated systematically by Kang [21-23]. His 
results showed that there are significant individual or combined 
effects of these factors on boiling heat transfer characteristics. 


2.2. Boiling heat transfer of nanofluids—literature survey 


The boiling critical or burnout heat flux enhancement of 
particle-suspensions depends on the particle concentration, pH of 
the solution as well as on the deposition of the particles on the sus- 
pended heater surface. It is long back proven that addition of solid 
particle in base fluid can alter its boiling heat transfer performance. 
For example, Yang and Maa [7] first used nano-sized Al203 particles 
as small as 50nm in water to study the pool boiling heat transfer 
characteristics. They reported significant increase in pool boiling 
performance for very small volumetric concentrations (0.1-0.5%) 
of nanoparticles. After nanofluids emerged a growing number of 
research groups have got involved with boiling heat transfer study 
of nanofluids and it is timely to review their research findings. 

For Al203-water nanofluid with a flat plate heater, You et al. 
[24] observed a three-fold increase in critical heat flux (CHF). For 
silica nanofluids, similar three-fold enhancement in CHF was later 
reported by Milanova and Kumar [25]. Whereas, Das et al. [26] 
reported deterioration of boiling heat transfer of water in the pres- 
ence of Al203 nanoparticles in it. Their outcome was partially 
attributed to the properties of the nanofluid, boiling surface and 
the interaction between the two. In contrast to Das et al. [26], Wen 
and Ding [27] showed that the enhancement of pool boiling heat 
transfer of the same y-Al)03/water-based nanofluid was about 40% 
at 1.25 wt.% of particle loading. In another study, Witharana [28] 
investigated the boiling heat transfer performance of two types 
of Au and SiOz-laden aqueous nanofluids in a cylindrical vessel 
under atmospheric pressure. They found that the boiling heat trans- 
fer increases and decreases for Au-nanofluid and SiO2-nanofluid, 
respectively. These conflicting results are not well-explained. 

Prakash et al. [29,30] performed two sets of experiments to 
quantify the effect of heater orientation and heater surface rough- 
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ness on pool boiling heat transfer of Alp03/water-based nanofluids. 
Their first set of experiment [29] with vertical tubular heaters of 
various surface roughness showed that while the rough heater 
surface increases heat transfer, smooth surface significantly dete- 
riorates the heat transfer. For example, for rough heater surface 
the heat transfer enhancement was about 70% at 0.5 wt.% concen- 
tration of alumina nanoparticles. Whereas, for smooth heater the 
heat flux reduction reaches up to 45% at a particle concentration 
of 2 wt.%. Their other set of experiment [30] was mainly to iden- 
tify the influence of heater orientation on the boiling heat transfer 
of the same nanofluid (Al,03/water). Their results demonstrated 
that horizontal and inclined heater orientations showed enhance- 
ment and deterioration of boiling heat transfer of this nanofluid, 
respectively. 

Soltani et al. [31] investigated the pool boiling heat trans- 
fer performance of Newtonian nanofluids under various heat flux 
densities. In their study, y-Al203 (20-30nm)/water and SnO2 
(55nm)/water-based nanofluids were used in a vertical cylin- 
drical glass vessel. Their results showed that except for low 
concentrations (>0.5 wt.%) of SnO2 nanoparticles, the boiling heat 
transfer coefficients of these nanofluids increase with increas- 
ing concentration of nanoparticles. These paradoxical results were 
attributed to the differences in thermal conductivity and size of 
these two nanoparticles. In a similar pool boiling experiment, this 
same research group [32] recently used non-Newtonian nanoflu- 
ids of y-Al,O3 nanoparticles dispersed in carboxy methyl cellulose 
(cmc)-water solution and the boiling heat transfer coefficient was 
found to increase by about 25% at 1.4 wt.% loading of nanoparticles. 
Recently, Truong et al. [33] conducted pool boiling experiments of 
diamond, ZnO, and Al203 nanoparticles-laden aqueous nanofluids 
with modification of sandblasted as well as bare plate heaters. They 
found up to 35% increase in CHF for pre-coated heaters compared 
to those of bare and sandblasted heaters. 

Among very few studies on flow boiling of nanofluids, Kim 
et al. [34] reported about 50% enhancement in flow boiling CHF for 
Al203/water nanofluids flowing through a vertical stainless steel 
tube. Very recently, Henderson et al. [35] studied refrigerant-based 
SiOz and CuO-nanofluids in flow boiling experiments in horizon- 
tal copper tube. They found that while the boiling heat transfer 
coefficient (BHTC) of SiO2/R-134a nanofluid decreases up to 55% in 
comparison to pure R-134a, the BHTC increases more than 100% 
for CuO-laden nanofluid over base fluid, i.e. mixture of R-134a and 
polyolester oil (PO). 

A summary of studies on boiling heat transfer of nanofluids is 
presented in Table 1. It can be noticed that besides inconsistent 
and contradictory results, most of the researchers used alumina- 
nanofluids. Whereas very few efforts have been made with carbon 
nanotubes (CNTs)-nanofluids which exhibit much higher thermal 
performance compared to those of ceramic-nanofluids [14,15]. Itis, 
therefore worthwhile to study boiling heat transfer characteristics 
of CNTs-nanofluids. 

A comparison of heat flux versus superheat results from var- 
ious groups is shown in Fig. 4. From these representative results 
(Fig. 4), it can clearly be seen that heat flux (also critical heat flux) 
data relative to superheat reported by various research groups vary 
widely. This is probably due to the differences in characterization of 
nanofluids, different size and concentration of nanoparticles used 
and different types of heaters used in various research groups. 
Despite of the fact that some researchers reported deterioration 
of boiling heat transfer of nanofluids, the significant increase in the 
critical heat flux of nanofluids is still undisputed. 


2.3. Boiling heat transfer—correlations 


Although there are many correlations [18,20] to predict 
the boiling heat transfer coefficient and critical heat flux, 


0.1 wt.% Al,O, in water [7] 
0.5 wt.% Al,O, in water [7] 
0.1 vol.% Al,O, in water [26] 
1 vol.% Al,O, in water [26] 
0.71 wt.% Al,O, in water [27] 
1.25 wt.% Al,O, in water [27] 
0.6 wt.% Al,O, in water [31] 
1 wt.% Al,O;, in water [31] 
0.5 vol.% Al,O, in water [37] 


1 vol.% Al,O, in water [37] 
0.15 vol.% ZrO, in water [41] 
0.1 wt.% CuO in water [43] 
0.5 wt.% CuO in water [43] 
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Fig. 4. Comparison of heat flux versus superheat results from various research 
groups. 


researchers working on boiling heat transfer with nanofluids 
[27,30,36,37] mainly used predictions from two most popular 
correlations—Rohsenow [45] and Zuber [46] correlations to com- 
pare their experimental data. Thus, without elaborating details of 
other models or correlations, these two classical correlations are 
presented here. 

Based on dimensional analysis of all relevant factors arising in 
nucleate boiling and experimental data over a wide variety of con- 
ditions Rohsenow [45] proposed one of the first and most widely 
used nucleate boiling correlations which is commonly expressed in 
the following form: 


P 0.33 n 
Cp(Tw — Ts) =c,|4 o (2) (1) 
hg Y| hg V 8(0 — bs) k 


where Cp, hfg, Tw, Ts, q”, H, P, Pg, 0, and k are specific heat of liquid, 
the latent heat of the fluid, temperature of heater wall, saturation 
temperature of liquid, heat flux, viscosity of saturated liquid, den- 
sity of saturated liquid, density of saturated vapor, surface tension, 
acceleration of gravity and thermal conductivity of liquid, respec- 
tively. The values for empirical coefficient C,; and exponent n for 
various surface-fluid combinations can be found elsewhere [47]. 
For example, Cyr and n for water-platinum combination are 0.013 
and 1, respectively. 

The maximum heat flux attainable in the nucleate pool boil- 
ing regimes is the critical heat flux (CHF) or point of departure 
from nucleate boiling (DNB). A widely used correlation which pre- 
dicts the pool boiling CHF for a flat infinite heater surface (facing 
upwards) is Zuber’s correlation [46]. The Zuber CHF is expressed 
as: 


diur = Kol’ hg [goo — pg)| (2) 

The value of constant K ranged from 0.138 to 0.157. However, 
simplifying the formulation analysis, Zuber proposed a value of 
K=0.131. Although widely used, most of the cases predictions of 
pool boiling heat transfer and critical heat flux of nanofluids using 
these two abovementioned correlations do not agree well with the 
experimental results [27,36]. There is therefore a strong need to 
develop new correlations or modify these classical models to be 
used for nanofluids. 


2.4. Present boiling experiments with CNT-nanofluids 
For pool boiling experiments, high purity single-walled carbon 


nanotubes (SWCNTs) are suspended in deionized water to prepare 
sample nanofluids. As a part of surface treatment, CNT bundles were 
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Summary of boiling experiments with nanofluids. 


Researchers 


Heater (type of boiling) 


Nanofluids 


Remarks 


Yang and Maa [7] 
Witharana [28] 
Das et al. [26] 

You et al. [24] 
Vassallo et al. [36] 
Bang and Chang [37] 
Wen and Ding [27] 
Kim et al. [38] 

Kim et al. [39] 
Jackson [40] 
Prakash et al. [29] 


Prakash et al. [30] 


Chopkar et al. [41] 


Liu and Liao [42] 


Lv and Liu [43] 


Kathiravan et al. [44] 


Kim et al. [34] 
Soltani et al. [31] 
Soltani et al. [32] 


Henderson et al. [35] 


Horizontal tube heater 


(pool) 
Cylindrical vessels (pool) 


Cylindrical cartridge (pool) 
Cartridge (pool) 

NiCr wire (pool) 

Square flat heater (pool) 
Flat disk heater (pool) 


Smooth NiCr wire heater 
(pool) 

Stainless steel wire and flat 
plate heaters (pool) 

Flat copper coupon heater 
(pool) 

Vertical tubular heaters 
(pool) 

Tubular heaters at various 
orientations (pool) 


Flat surface (pool) 


Flat horizontal copper 
surface (pool) 


Vertical small heated tubes 
(pool) 


Horizontal tube (pool) 


Vertical stainless steel tube 
(flow) 

Vertical cylindrical glass 
vessel (pool) 

Vertical cylindrical glass 
vessel (pool) 

Horizontal copper tube 
(flow) 


Al203/water 

Au and SiO2/water and EG 
Al203/water 

Al203/water 

SiO2/water 

Al203/water 

Al203/water 

TiO2/water 

Al203, ZrOz2 and SiO2/water 
Au/water 


Al203/water 


Al203/water 


ZrO2/water 


CuO/water and Alcohol 
(C2Hs50H) 


CuO/water 


CNTs/water 


Al203/water 
Al203 and SnO2/water 
Al203/cmc-water solution 


SiO2/R-134a and 
CuO/R134a + polyolester oil 


Boiling heat flux increases considerably. 


Boiling heat transfer increases for Au-nanofluids but 
decreases for SiO2-nanofluids. 

Deterioration of heat transfer performance. 

Critical heat flux (CHF) increases up to 200%. 

CHF increases significantly. 

Pool boiling heat transfer deteriorates but CHF increases. 
The boiling heat transfer increases up to 40% at 1.2 wt.% of 
nanoparticles. 

The maximum enhancement of CHF is 200%. 


The maximum increase in CHF is 80% but the heat transfer 
rate is deteriorated. 

While the heat transfer decreased about 20%, the 
maximum CHF increase is 5 times over water. 

While rough heater surface increases heat transfer, smooth 
surface significantly deteriorates. 

Horizontal and inclined heater orientations showed 
enhancement and deterioration of heat transfer, 
respectively. 

Enhanced boiling heat transfer is found at low particle 
concentration. Tetramethyl ammonium hydroxide 
surfactant increases the heat transfer rate. 

The boiling heat transfer characteristics of the both 
nanofluids are somewhat poor and the CHF values are 
higher than those of the base liquids. 

Enhancement in CHF and boiling heat transfer rate are 
observed for surfactant-free nanofluids but addition 
NaDBS surfactant deteriorates the heat transfer. 

Boiling heat transfer coefficient (BHTC) increases up to 
1.75-fold for 0.25 vol.% of CNT. Sodium lauryl sulfate (SLS) 
surfactant of 9 wt.% was also added. 

Flow boiling CHF increases about 50%. 


Except for low concentrations (>0.5 wt.%) of SnO2, the 
BHTC increases with loading of nanoparticles. 

The BHTC increases by 25% for 1.4 wt.% of Al203 
nanoparticles. 

While BHTC of Si02/R-134a nanofluid decreases up to 55% 
in comparison to pure R-134a, it increases more than 100% 


Sandblasted and bare 
horizontal plate heaters 
(pool) 


Truong et al. [33] 
Al203/water 


Diamond, ZnO and 


for CuO-laden nanofluid over base fluid. 
The CHF of the pre-coated heaters increased by up to 35% 
with respect to that of the bare, sandblasted heaters. 


refluxed with hydrochloric acid at 100 °C for several hours. This acid 
was chosen as a reactive reagent because it removes catalytic par- 
ticles without reducing the length of the tubes or damaging the 
side walls. Other acids (HNO3 and HNO3:H2SO,) tried as refluxing 
mixtures, were too harsh and damaged the structure. CNTs were 
separated from the acid by centrifugation for about an hour and 
left to dry. Different concentrations of sodium dodecyl benzene- 
sulfonate (NaDBS) surfactant are used as dispersing and stabilizing 
agent for nanotubes in water. The schematic of experimental facili- 
ties used is shown in Fig. 5. Details of the experimental facilities and 
procedures are reported in previous studies [48,49] and will not be 
discussed here. Representative results of surfactant- and surface 
tension-dependent pool boiling experiments of CNTs-nanofluids 
[48] are presented and discussed in this section. Both the critical 
heat flux (CHF), which is the sudden jump in temperature at one 
and the same heat flux and the burnout heat flux (BHF), which is the 
case when the wire is brought to catastrophic failure (i.e. complete 
burnt out) are determined at constant 0.1% volumetric loading of 
CNT and for various concentrations of surfactant. 


2.4.1. Effect of surfactant concentration on boiling heat flux 

The results of pool boiling heat flux with respect to superheat 
(wire temperature minus saturation temperature of liquid, Tw — Ts) 
are presented and compared with the one for pure de-ionized water 


in Fig. 6. The NaDBS surfactant to CNTs concentrations are varied 
from 1:20 to 1:1. Fig. 6 demonstrates that the CNT-nanofluid with 
any concentration of surfactant exhibits higher CHF value than that 
of base fluid. The CHF value of pure deionized water is 750 kW/m2. 
The effect of increasing the surfactant concentration from 1:20 to 
1:5 (NaDBS:CNT) is that the CHF increases as well. However, if the 
concentration of surfactant is further increased from 1:5 to 1:1, the 
CHF drops drastically from 4439 kW/m? to 1322 kW/m2, respec- 
tively. There is, therefore a critical concentration of surfactant for 
which the CHF and BHF reach maximum values. The highest CHF 
and BHF results are obtained for the concentration ratio of 1:5. 
The deposition of nanoparticles on the heater wire is believed to 
be one of the key reasons for any enhancement in CHF of this 
nanofluid up to critical concentration of surfactant. Kim et al. [50] 
also claimed that the deposition of nanoparticles on the heater wire 
is the main reason for the enhancement of CHF of their nanofluids. 
SEM micrographs shown in Fig. 7 demonstrate that with an increase 
in surfactant the deposition on the wire increases as well. However, 
the maximum CHF and BHF do not occur at the maximum particle 
deposition, i.e. at 1:1 NaDBS:CNT ratio. 


2.4.2. Effect of surface tension on burnout heat flux 
The burnout heat flux is determined at constant 0.1% volumetric 
loading of CNT and at various surface tensions of nanofluids which 


S.M.S. Murshed et al. / Renewable and Sustainable Energy Reviews 15 (2011) 2342-2354 2347 


ARP 
pem 


Water in ——4 ħ—>Water out 


Coiled 
condenser 


Multimeter 
Power supply 


SE 


PC with 
DAQ 


Ñ 


Lead wires 


Nanofluid 
Thermocouple 


Ni-Cr wire 


—— Heater 


Fig. 5. Schematic of pool boiling experimental setup [48]. 
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Fig. 6. Effect of surfactant concentration on heat flux of SWCNTs-nanofluid. 


was adjusted by varying concentrations of NaDBS surfactant. The 
surface tensions of nanofluid with surfactant (oys) and deionized 
water with surfactant (ows) concentrations are measured. At any 
concentration of surfactant, the difference between surface ten- 
sions of nanofluid and its based fluid (deionized water) is termed 
as relaxation of surface tension, i.e. Ao; = Oys — ows. This relaxation 
of surface tension of nanofluids is the driving force for the burnout 
heat flux of nanofluid, i.e. (BHF)y in pool boiling as shown in Fig. 8. 
Based on the best fit of the experimental data, a nonlinear empir- 
ical correlation between the BHF and relaxation of surface tension 
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Fig. 8. Effect of surface tension on the burnout heat flux of nanofluid (NaDBS to 
CNTs concentration ratios are shown in parentheses). 


is obtained and has the form 


(BHF) _ 4 + 0.073 A0; — 0.0065Ao;? + 0.0003 A07? (3) 
(BHF)w 

where (BHF)w is the burnout heat flux of water and 
Ao; =(Ons — ows) in mN/m. 

It is demonstrated that regardless of surfactant concentration, 
the burnout heat flux of nanofluid strongly depends (i.e. nonlinearly 
related) on its relaxation of surface tension with base fluid. Similar 
to the CHF, the BHF of nanofluid also reaches its highest value at 
NaDBS and CNT concentration ratio of 1:5. The highest BHF value of 
nanofluid is about 265% enhancement over that for the base fluid at 
Aor 25 mN/m. The burnout heat flux of deionized water (BHF)w is 
found to be 1500 kW/m2. The observed increment of relative BHF of 
nanofluid (i.e. (BHF)y/(BHF)w) can be attributed to the Marangoni 
effect which results from the surface tension differences. Both the 
CHF and burnout are found to happen at the optimum 1:5 sur- 
factant:CNT ratio. Thus, the pool boiling heat transfer behavior of 
surfactant-added nanofluid is not only dictated by the deposition 
of nanotubes on heater surface, but also by the relaxation of sur- 
face tension of nanofluid, which is a precursor to the deposition and 
overall onsets of both the CHF and BHF. 


2.5. Studies on nanofluids droplet spreading 


Droplets impingement and spreading studies on heated sub- 
strate surface are of great importance for the practical application 
of nanofluids as an advanced coolant particularly in spray cooling. 
Unfortunately, very little research efforts have been made to study 
nanofluids’ droplet impingement behavior on solid surfaces under 
various conditions. Among very few studies, Wasan and Nikolov 
[51] were the first to investigate the effects of the particle structure 
formation and the structural disjoining pressure of nanoparti- 
cles on the spreading of nanofluids on solid surface. Recently, 
Duursma et al. [52] studied the effect of aluminum nanoparti- 
cles on droplets boil-off by allowing nanofluid drops to fall onto 
a copper surface at temperature higher than the liquid satura- 
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Fig. 7. SEM images of NiCr wires at various surfactant concentrations after pool boiling experiments [48]. 


2348 S.M.S. Murshed et al. / Renewable and Sustainable Energy Reviews 15 (2011) 2342-2354 


tion temperature. They demonstrated that increasing the surface 
temperature and Weber number promote the receding breakup 
scenario, while increasing the nanoparticle concentration discour- 
ages this breakup. In another recent study, the influence of surface 
temperature on the hydrodynamic characteristics of water and a 
nanofluid droplets impinging on a polished and nanostructured 
surface was investigated by Shen et al. [53]. Their results showed 
that SWCNTs-nanofluid has larger spreading diameter compared to 
that of deionized water and using a nanofluid or a nanostructured 
surface can reduce the total evaporation time up to 37%. Murshed 
and Nieto de Castro [54] very recently reported the spreading char- 
acteristics of ananofluid droplets impinging on a metallic substrate 
under the influence of several key factors such as nanoparticles 
volume fraction, substrate temperature, and the Weber number. 
Results showed that the transient spreading diameter and height 
of droplet impacting onto solid surface are greatly influenced by 
each of these factors. Nevertheless, more studies are needed on 
dynamics of non-boiling and boiling nanofluid droplets impinging 
on solid or liquid surfaces as the spreading of liquid droplet plays 
a key role in many industrial processes like spray cooling, coating, 
ink-jet printing, and oily soil removal. 


3. Flow and heat transfer characteristics of nanofluids 


Compared to the reported studies on thermal conductiv- 
ity, investigations on convective heat transfer of nanofluids are 
still scarce. However, the practical applications of nanofluids as 
advanced heat transfer fluids are mainly in flowing systems such as 
mini- or micro-channels heat sinks and miniaturized heat exchang- 
ers. In this section, we critically review both experimental and 
computational studies on convective heat transfer of nanofluids. 
In addition, some representative results from our experimental 
investigation on laminar flow convective heat transfer of TiO2- 
nanofluids are also presented and discussed. 


3.1. Experimental studies on convective heat transfer of 
nanofluids 


The first experiment on convective heat transfer of nanofluids 
(e.g. y-Al203/water) under turbulent flow conditions was per- 
formed by Pak and Cho [55]. In their study, even though the Nusselt 
number (Nu) was found to increase with increasing nanoparticle 
volume fraction and Reynolds number, the heat transfer coefficient 
(h) actually decreased by 3-12%. The reasons for such paradoxical 
results might be the observed large enhancement in viscosity. On 
the other hand, Eastman et al. [56] later showed that with less than 
1vol.% of CuO nanoparticles, the convective heat transfer coeffi- 
cient (h) of water increased more than 15%. The results of Xuan 
and Li [57] illustrated that the Nusselt number of Cu/water-based 
nanofluids increased significantly with the volumetric loading 
of particles as for 2vol.% of nanoparticles, the Nusselt number 
increased by about 60%. Wen and Ding [58] investigated the heat 
transfer behavior of nanofluids at the tube entrance region under 
laminar flow conditions and showed that the local heat transfer 
coefficient varied with particle volume fraction (@) and Reynolds 
number (Re). They also observed that the enhancement is par- 
ticularly significant at the entrance region. Later Heris et al. [59] 
studied convective heat transfer of CuO and Al203/water-based 
nanofluids under laminar flow conditions through an annular tube. 
Their results showed that heat transfer coefficient increases with 
particle volume fraction as well as Peclet number. In their study, 
Al,O3/water-based nanofluids found to have higher enhancement 
of heat transfer coefficient compared to CuO/water-based nanoflu- 
ids. 
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Fig. 9. Comparison convective heat transfer results from various research groups 
[55,57,58,62,70]. 


An experimental investigation on the forced convective heat 
transfer and flow characteristics of TiO -water nanofluids under 
turbulent flow conditions is reported by Duangthongsuk and 
Wongwises [60]. A horizontal double-tube counter flow heat 
exchanger is used in their study. They observed a slightly higher 
(6-11%) heat transfer coefficient for nanofluid compared to pure 
water. The heat transfer coefficient increases with increasing mass 
flow rate of the hot water as well as nanofluid. They also claimed 
that the use of the nanofluid has a little penalty in pressure drop. 

In microchannel flow of nanofluids, Faulkner et al. [61] was the 
first to conduct convective heat transfer experiments with aque- 
ous CNT-nanofluid in a microchannel with hydraulic diameter of 
355 um at Reynolds numbers between 2 and 17. They found con- 
siderable enhancement in heat transfer coefficient of this nanofluid 
at CNT concentration of 4.4%. Later, Jung et al. [62] studied heat 
transfer performance of Alj03/water-based nanofluid in a rectan- 
gular microchannel under laminar flow condition. Results showed 
that the heat transfer coefficient increased by more than 32% for 
1.8 vol.% of nanoparticles and the Nusselt number increases with 
increasing Reynolds number in the flow regime of 5 > Re < 300. 

An up-to-date overview of the published experimental results 
on the convective heat transfer characteristics of nanofluids is 
given in Table 2. A comparison of results of Nusselt number versus 
Reynolds number for both laminar and turbulent flow conditions 
from various groups is shown in Fig. 9. From Table 2 and Fig. 9, it 
can be seen that the results from various groups vary widely and 
most of the studies lack physics-based explanation of the observed 
results. Although there are couple of review articles [71,72] that 
attempted to put together some studies on convective heat trans- 
fer of nanofluids, no comprehensive analysis of up-to-date findings 
are reported. 

Among very few research efforts, the work of Putra et al. [73] 
was the first to investigate the natural convective heat transfer of 
aqueous CuO and Al,03-nanofluids. They used a horizontal poly- 
oxymethylene cylinder which was heated from one end and cooled 
from the other. Significant deterioration (decrease) of convection 
heat transfer for these nanofluids was observed and the deteriora- 
tion increases with increasing particle concentration, particularly 
for CuO nanoparticles. Putra et al. ascribed the possible reasons of 
such deterioration to the effects of particle—fluid slip and sedimen- 
tation of nanoparticles. In contrast to Putra et al. [73], numerical 
simulation of Khanafer et al. [74] showed that in a 2D horizon- 
tal enclosure, the natural convection heat transfer coefficient (Nu) 
of nanofluids increases with particle concentration. By taking into 
account the solid particle dispersion, they [74] developed a model 
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Table 2 
Summary of forced convection heat transfer experimental studies of nanofluids. 
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Researchers Geometry/flow nature Nanofluids 


Findings 


Pak and Cho [55] Tube/turbulent 


Xuan and Li [57] Tube/turbulent Cu/water 
Wen and Ding [58] Tube/laminar Al203/water 
Ding et al. [63] Tube/laminar CNT/water 


Yang et al. [64] Tube/laminar 


Heris et al. [59] Tube/laminar 


Lai et al. [65] 
Jung et al. [62] 


Tube/laminar 
Rectangular 
microchannel/laminar 
Tube/turbulent 
Tube/laminar 


Williams et al. [66] 
Hwang et al. [67] 


Xie et al. [68] Tube/laminar Al203, ZnO, TiO and 
MgO/water 
Amrollahi et al. [69] Tube/laminar and MWCNT/water 


turbulent 


Al203 and TiO2/water 


Graphite/automatic 
transmission fluid 
Al203 and CuO/water 


Al203/water 
Al203/water 


Al203 and ZrO2/water 
Al203/water 


At 3 vol %, h was 12% smaller than pure water for a given 
average fluid velocity. 

A larger enhancement of h with volume fraction (ġ) and Re 
was observed. 

Increased h with ġ and Reynolds number was observed. 
NaDBS was used as the surfactant 

At 0.5 wt.%, h increased by more than 350% at Re = 800. 
The nanoparticles considerably increase the heat transfer 
coefficient of the fluid system in laminar flow. 

h increase with @ and Pe. Al203 shows higher 
enhancement than that of CuO. 

Nu increased 8% for @=0.01 and Re=270. 

h increased 15% for #=0.018. 


h increased significantly. 
h increased only up to 8% at Re =730 for #=0.003. 
h increased up to 252% at Re = 1000 for MgO nanofluid. 


h increased up to 33-40% at concentration of 0.25 wt.%. 


to analyze the heat transfer performance of nanofluids inside an 
enclosure. 

Wen and Ding [75] later investigated heat transfer behavior of 
specially formulated TiOz/water nanofluid under the natural con- 
vection conditions. The results showed that the natural convective 
heat transfer coefficient decreases with increasing particle con- 
centration. These unexpected results are in contradiction to the 
numerical findings of Khanafer et al. [74] but are in agreement with 
the observations by Putra et al. [73]. Along with the reasons ascribed 
by Putra et al. [73], Wen and Ding [75] added several other rea- 
sons for such paradoxical results that include convection induced 
by concentration difference, modification of dispersion properties 
and particle-surface and particle-particle interactions. 


3.2. Theoretical models for convective heat transfer of nanofluids 


Researchers investigating convective heat transfer of nanoflu- 
ids employed existing conventional single-phase fluid correlation 
to predict the heat transfer coefficient or proposed new correlations 
obtained by fitting their own experimental data. However, none of 
these correlations were validated with wide range of experimental 
data under various conditions and thus are not widely accepted. In 
an attempt to establish a strong explanation of the reported anoma- 
lously enhanced convective heat transfer coefficient of nanofluid, 
Buongiorno [76] considered seven-slip mechanisms and concluded 
that among those seven only Brownian diffusion and thermophore- 
sis are the two most important particle/fluid slip mechanisms in 
nanofluids. Besides proposing a new correlation, he also claimed 
that the enhanced laminar flow convective heat transfer can be 
attributed to a reduction of viscosity within and consequent thin- 
ning of the laminar sublayer. Commonly used classical models 
together with recently proposed correlations for the predictions 
of convective heat transfer coefficient of nanofluids are summa- 
rized in Table 3. Considering all possible key mechanisms such as 
particle Brownian motion and migration, thermophoresis and iner- 
tial a generalized convective heat transfer correlation needs to be 
developed. 


3.3. Numerical studies on convective heat transfer of nanofluids 


Several preliminary numerical studies on the convective heat 
transport of nanofluids were conducted by Khanafer et al. [74] 
and a research group from Université de Moncton, Canada [79,80]. 
Results from the Canadian group showed that nanofluids can sub- 


stantially enhance heat removal in both tube [79] and radial flows 
[80]. For example, Roy et al. [80] numerically studied the hydro- 
dynamic and thermal fields of Alz03/water-based nanofluids in a 
radial laminar flow system and reported a two-fold increase in heat 
transfer for 10 vol.% of nanoparticles in water. They also observed 
that the wall shear stress increased with an increase in particle vol- 
ume fraction. For Aly03/water- and Al203/ethylene glycol-based 
nanofluids, numerical results of Maiga et al. [79] showed that the 
heat transfer coefficient increases considerably with increasing 
nanoparticle concentration. However, the presence of such parti- 
cles has also induced drastic effects on the wall shear stress that 
increases appreciably with the particle loadings. 

A Lattice Boltzmann model for simulating flow and energy trans- 
port processes inside nanofluids was proposed by Xuan et al. [81]. 
They considered the external and internal forces acting on the sus- 
pended nanoparticles and interactions among the nanoparticles 
and fluid particles. Based on the simulation results, they inferred 
that the fluctuation of Nusselt number of nanofluids along the main 
flow direction is due to the unstable distribution of nanoparticles. 
No further development is made in exploring the use of the Lattice 
Boltzmann method to explain the anomalous thermal convection 
of nanofluids. 

Palm et al. [82] performed numerical simulations for the laminar 
forced convection flow of nanofluids with temperature-dependent 
properties. They found that Al,03/water-based nanofluids with 
a volume fraction of 4% can produce a 25% increase in the 
average wall heat transfer coefficient compared to water alone. 
Significant differences were also found when using constant 
property nanofluids (temperature-independent) versus nanofluids 
with temperature-dependent properties. The use of temperature- 
dependent properties resulted in larger heat transfer predictions 
with corresponding decrease in wall shear stresses when com- 
pared to predictions using constant properties nanofluids. With an 
increase in wall heat flux, they found that the average heat transfer 
coefficient increases while the wall shear stress decreases. 

Based on Khanafer et al.’s [74] model and taking into account 
particle dispersion, Jou and Tzeng [83] numerically studied 
convective heat transfer performance of nanofluids inside two- 
dimensional rectangular enclosures. Their results showed that 
increasing volume fraction and buoyancy parameter can cause an 
increase in average heat transfer coefficient. 

The cooling performance of a microchannel (silicon) heat sink 
under forced convective laminar flow with nanofluids is numeri- 
cally investigated by Jang and Choi [84]. Results showed that at a 
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Table 3 


Correlations used for predicting the convective heat transfer coefficient of nanofluids. 


Origin/reference Correlations 


Flow regimes/conditions/remarks 


Shah equation [77] Nu = 


1/3 
1.953 (ser?) : 


4.364 + 0.0722Re Pr x i 


Dittus and Boelter [78] Nu=0.023Re°8Ppr°4 


Sieder-Tate equation [77] Nu =0.027Re°8 Pr" (j2,/ us )°'4 


where up and us are the fluid viscosity at bulk 


for (rer?) > 33.3 


for RePr z < 33.3 


Laminar flow and constant heat flux condition. It is popular 
for a thermal entrance region. 


Heating of fluid and used for 0.7 < Pr < 100. It is very 
popular in pipe flow cases and is not recommended for 
temperature dependent property variation. 

Valid for 0.7 < Pr < 16,700 and Re > 10,000. It is used for 
temperature-dependent property variation cases. 


fluid temperature and at heat-transfer boundary 


surface temperature, respectively. 
Nu=0.021Re®8Ppr°5 
Laminar flow: 


Pak and Cho [55] 
Xuan and Li [57] 


Nu =0.4328(1 + 11.2850754 Pe®?18 )Re0-33 p04 


Turbulent flow: 


Turbulent flow 
Pe is the particle Peclet number. It considered particles 
concentration. 


Nu =0.0059(1 + 7.6286¢°-5886 pe0.001 )Re0.9238 py04 


Maiga et al. [79] Averaged Nu for constant heat flux: 

Nu = 0.086Re°* pr’? 

Averaged Nu for constant temperature: 
Nu = 0.28Re°® pr°6 


Jung et al. [62] Nu=0,0095p°! Re®-42 p55 


(f/8)(Rep —1000)Pry, 


Buongiorno [76] Nu = ot Siar 1 


Correlations are proposed from numerical study and valid 
for Re < 1000 and 6 < Pr < 753. 


For microchannel laminar flow regime. It considered the 
volume fraction of nanoparticles. 


Turbulent flow 


where f is the friction factor and ôf is an empirical 


constant. 


fixed pumping power (2.245 W) the cooling performance of their 
heat sink with aqueous-diamond (1 vol.%) nanofluid is enhanced 
by about 10% compared to that of with water alone. They also 
reported that nanofluids reduced both the thermal resistance and 
the temperature difference between the heated microchannel wall 
and coolant. 

Using a dispersion model Heris et al. [85] later performed 
numerical simulations for the laminar-flow convective heat trans- 
fer of aqueous nanofluids in a circular tube with constant wall 
temperature. Effects of concentration and size of various types 
of nanoparticles (Al,03, CuO, and Cu in water) as well as Peclet 
number on heat transfer coefficient are investigated. They reported 
that Nu increases with increasing particle concentration and with 
decreasing nanoparticles size. Their numerical predictions were 
found to be in considerable agreement with their experimental 
data. 

A numerical study on laminar mixed convection of Alz03/water 
nanofluid was reported by Mirmasoumi and Behzadmehr [86]. 
Two-phase mixture model was used to study the effects of nanopar- 
ticles mean diameter on the flow parameters. Their calculated 
results showed that like Heris et al.’s [85] findings, the heat transfer 
coefficient significantly increases with decreasing the nanoparti- 
cles mean diameter. However particle size does not notably change 
the skin friction coefficient. 

Ina separate numerical study, Ho et al. [87] attempted to identify 
the effects of uncertainties in various effective viscosity and ther- 
mal conductivity models adopted on buoyancy-driven convection 
heat transfer of alumina—water nanofluid in a square enclosure. 
Based on their simulation results from a comparative study of four 
different models, they demonstrated that the effect of nanopar- 
ticle concentration on the averaged Nusselt number significantly 
depends on the effective viscosity and thermal conductivity models 
used. 

Oztop and Abu-Nada [88] recently numerically investigated 
the buoyancy-driven natural convection heat transfer character- 
istics of different types of water-based nanofluids containing Cu, 
Al,03, and TiO nanoparticles in a partially heated enclosure. An 
increase in mean Nusselt number was observed with the volume 


fraction of nanoparticles for the whole range of Rayleigh number 
(103 <Ra<5 x 10°) and heat transfer also increases with increasing 
of height of heater. 

A very recent numerical work on a 3D laminar flow heat trans- 
fer with Al203, CuO/mixture of ethylene glycol and water-based 
nanofluids circulating through flat tubes of a radiator was under- 
taken by Vajjha et al. [89]. Their results demonstrated that the 
convective heat transfer coefficient in the developing and devel- 
oped regions along the flat tubes with the nanofluid flow showed 
remarkable improvement over the base fluid and for the local 
and the average friction factor and convective heat transfer coef- 
ficient show an increase with increasing volumetric concentration 
of nanoparticles. 

Above review on numerical studies demonstrated that convec- 
tive heat transfer coefficient of nanofluids considerably increases 
with increasing concentration of nanoparticles and with decreas- 
ing nanoparticle size. However, more works on forced and natural 
convective heat transfer of nanofluids are to be performed by using 
advanced numerical methods such as Lattice Boltzmann, Molecular 
and Brownian dynamic simulations. 


3.4, Present experiments on laminar flow convective heat transfer 
of nanofluids 


Sample nanofluids for convective heat transfer experiments 
were prepared by dispersing different volume percentages, i.e. 
0.2-0.8% of titanium dioxide (TiO ) nanoparticles (15 nm diameter) 
in deionized water (DIW). As a dispersant agent, Cetyl Trimethyl 
Ammonium Bromide (CTAB) surfactant was added in the sample. 
To ensure proper dispersion of nanoparticles, sample nanofluid was 
homogenized by using an ultrasonic dismembrator. 

The effects of nanoparticles concentration and Reynolds number 
on the convective heat transfer coefficient of TiO2/deionized water 
based-nanofluids were previously investigated [90] and some rep- 
resentative results are presented. An experimental setup was 
established to conduct experiments on heat transfer of nanofluids 
at laminar flow regime in a cylindrical channel. The schematic of 
experimental facilities used is shown in Fig. 10. Details of the exper- 
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Fig. 10. Schematic of convective heat transfer experimental setup [90]. 


imental facilities and procedures are reported elsewhere [90] and 
will not be discussed here. Instead formulations for obtaining the 
heat transfer coefficient (h) together with a most commonly used 
existing correlation for predicting the Nusselt numbers of nanoflu- 
ids in laminar flow conditions are presented. 

As detailed in previous paper [90], applying first law (energy bal- 
ance) in control volume of flow channel, the following formulation 
for the local heat transfer coefficient is obtained 


Pap. = 


q' 
(To,w(x) — ((q[2D3 In(Do/D;) — (D3 — D? )])/47(D4 — D? )ksx)} — (T; + ((To — T;)/L)x} 


(4) 


where To w(x) is the outer wall temperature of the tube (mea- 
surable), q” is the heat flux of the test section (W/m2), q is the heat 
supplied to the test section (W), ks is the thermal conductivity of the 
copper tube (W/m K), D; and Do are the inner and outer diameters 
of the tube, respectively and x represents the longitudinal location 
of the section of interest from the entrance. L is the length of the 
test section and T; and To are the inlet and outlet fluid temperature, 
respectively. 

Once the local heat transfer coefficient is determined and the 
thermal conductivity of the medium is known, the local Nusselt 
number is calculated from 


hnf-xDi 


= (5) 


Nung-x = 


where kyy is the effective thermal conductivity of nanofluids. The 
classical Hamilton-Crosser model is used for the determination of 
knf which is given by [91] 


kp + (n — 1)kp — (n 
knf = kf 7 kp i (n f 


1)Ø(kf — kp) 
1)ky t Pky kp) 


where ky and kp are the thermal conductivities of the base liquid and 
the nanoparticles, respectively, ¢ is the volume fraction of nanopar- 
ticles and n is the empirical shape factor, which has a value of 3 for 
spherical particle. 

The Nusselt number can also be determined from the existing 
correlations. The well-known Shah’s correlation for laminar flows 
under the constant heat flux boundary conditions is used and repro- 
duced as [77] 
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Fig. 11. Axial profiles of local heat transfer coefficient of nanofluid at Re = 1100 [90]. 


For steady and incompressible flow of nanofluids in a tube of 
uniform cross-sectional area, the Reynolds number and Prandtl 
number are defined as follows 


ail and Pr= onf Hf 


Re= 
TDi Lng knf 


(8) 


where m is the mass flow rate and [yf and Cp-nf are the viscosity 
and specific heat of nanofluids, respectively. 

While the specific heat of nanofluids is calculated using the fol- 
lowing volume fractioned-based mixture rule [55,79,84] 


Cy-nf = PCp-p + (1 m= P)Cp-f> (9) 


the viscosity of nanofluids is determined from Batchelor’s model 
given by [92] 


Unf = Uf + 2.56 + 6.247) (10) 


where @ is particle volume fraction and upis the base fluid viscosity. 
It is noted that other classical models for calculating the viscosity 
of mixture also yield similar results [13]. 


3.4.1. Axial profiles of local heat transfer coefficient 

Fig. 11 illustrates the local heat transfer coefficient against the 
axial distance from the entrance of the test section at Reynolds 
number (Re) of 1100. The results show that nanofluids exhibits con- 
siderably enhanced convective heat transfer coefficient which also 
increases with volumetric loadings of TiO2 nanoparticles. For exam- 
ple, at 0.8 vol.% of nanoparticles and at position x/D;=25 (where 
tube diameter D; = 4mm), the local heat transfer coefficient of this 
nanofluid was found to be about 12% higher compared to deionized 
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Fig. 12. Comparison with Shah’s correlation along axial distance at Re=1100 and 
1700 [90]. 
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Fig. 13. Effect of Reynolds number and nanoparticle volume fraction on Nusselt 
number at axial location of x/D; =25 [90]. 


water at the same Re. The enhancement in heat transfer coefficients 
of nanofluids is believed due to the enhanced effective thermal con- 
ductivity and the acceleration of the energy exchange process in the 
fluid due to the random movements of the nanoparticles. Another 
reason for such enhancement can be the migration of nanoparticles 
in base fluids due to shear action, viscosity gradient and Brownian 
motion in the cross section of the tube. 

Fig. 12 compares experimentally determined Nusselt numbers 
with the predictions by Shah’s correlation, i.e. Eq. (7) along the 
axial distance. Shah’s correlation slightly over-predicts the present 
results. The difference in tube size may be one of the reasons for 
such over prediction. A relatively small tube (4mm diameter) was 
used in this experiment, whereas the Shah’s equation was devel- 
oped on the basis of laminar flow in large channels. Nevertheless, 
similar over prediction by Shah’s equation was also reported by 
Wen and Ding [58]. 


3.4.2. Effect of Reynolds number and particle volume fraction on 
Nusselt number 

The effect of Reynolds number on Nusselt number is shown 
in Fig. 13. It can be seen that the measured Nusselt numbers for 
nanofluids are higher than those of water and they increase remark- 
ably with Reynolds number. The observed enhancement of the 
Nusselt number could be due to the suppression of the boundary 
layer, viscosity of nanofluids as well as dispersion of the nanoparti- 
cles. Fig. 13 also demonstrates particle volume fraction dependence 
of Nusselt number. The Nusselt number of this nanofluid is found 
to increase almost linearly with the particle volume fraction. The 


nanofluid behaves more like a fluid than a conventional solid 
(micrometer or millimeter)-fluid mixture. The effects of several 
factors such as gravity, Brownian force, and friction force between 
the fluid and the ultra-fine particles may coexist in the main flow 
of nanofluids. 


4. Concluding remarks 


A comprehensive review on three major cooling 
features—boiling, droplet spreading and convective heat transfer 
of nanofluids together with representative results from own 
experimental investigations on these areas are presented and 
analyzed in this paper. 

From the review of available experimental results for boiling 
heat transfer, it can be conferred that despite of contradictory 
and inconsistent data on boiling heat transfer, there is undisputed 
substantial increase in the boiling critical heat flux of nanofluids. 
However, reported data are still limited and scattered to clearly 
understand the underlying mechanisms and trend of boiling heat 
transfer performance of nanofluids. In addition, only a couple of 
efforts are made on flow boiling of nanofluids. Thus, it is impera- 
tive to conduct more research on flow boiling of nanofluids under 
the influence of various factors such as pressure, mass flux and 
subcooling besides performing more systematic experimental and 
theoretical investigations on pool boiling features of nanofluids. 
Without providing physical-chemistry based details explanation, 
there are common presumptions such as deposition of nano- 
particles or tubes on heat transfer surface and surface wettability 
used as reasons for the observed results on boiling heat transfer of 
nanofluids. 

This paper also presents investigations of the effects of surfac- 
tant concentration and surface tension on the pool boiling heat 
transfer of carbon nanotubes-nanofluid and results showed that 
large enhancement of boiling heat flux is possible and would 
depend on the concentration of the surfactants. The effect of relax- 
ation of surface tension of nanofluid on its burnout heat flux is found 
to be significant. 

Although spreading of liquid droplet plays a key role in many 
industrial processes like spray cooling, coating and ink-jet printing, 
only a couple of research efforts have been made to study droplet 
impingement characteristics of nanofluids on solid surfaces under 
various conditions. Thus more studies are needed on dynamics of 
non-boiling and boiling nanofluids droplets impinging on solid or 
liquid surfaces. 

The representative results and the review of the findings 
from the literature on convective heat transfer of nanofluids 
demonstrated that nanofluids exhibit an enhanced heat transfer 
coefficient compared to its base fluid and it increases significantly 
with increasing concentration of nanoparticles as well as Reynolds 
number. The review also shows a considerable chaos and random- 
ness in the reported data on convective heat transfer coefficient 
from various research groups. Therefore, more careful and system- 
atic investigations on nanofluids preparation and the convective 
heat transfer measurements are needed. A clear understanding of 
the convective heat transfer mechanisms in nanofluids is also not 
yet established. 

Researchers working on nanofluids have shown tremendous 
attention to the static thermal properties, particularly thermal 
conductivity of nanofluids rather than the boiling and convective 
heat transfer features. Although the applications of nanofluids as 
advanced cooling medium appears to be promising, the advance- 
ment toward concrete understanding on observed properties and 
features of nanofluids as well as their development for commercial 
applications remain challenging mainly due to lack of agreement 
in the data from different research groups, lack of understanding 
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of the mechanisms, and unsystematic measurements and sample 
preparation. Hence, proper sample preparation and repeatable and 
more systematic experimental studies on measuring any properties 
of nanofluids are worthwhile. 

In renewable energy sector, nanofluids could potentially be used 
to enhance heat energy storage from solar collectors and to increase 
the energy density. Thus, future research should also extend to such 
renewable energy-based applications of nanofluids. 


References 


[1] Lide DR. CRC handbook of chemistry and physics. Taylor and Francis; 2007. 

[2] Kaviany M. Principles of heat transfer. John Wiley & Sons, Inc.; 2002. 

[3] Masuda H, Ebata A, Teramae K, Hishinuma N. Alteration of thermal conductivity 
and viscosity of liquid by dispersing ultra-fine particles (dispersion of y-Al203, 
SiO2, and TiO» ultra-fine particles). Netsu Bussei 1993;4:227-33. 

[4] Slack GA. Thermal conductivity of MgO, AlzO3, MgAl2O4, and Fe30,4 Crystals 
from 3 to 300K. Physical Review 1962;126:427-41. 

[5] Slack GA. Thermal conductivity of pure and impure silicon, silicon carbide and 
diamond. Journal of Applied Physics 1964;35:3460-6. 

[6] Kim P, Shi L, Majumdar A, McEuen PL. Thermal transport measurements of 
individual multiwalled nanotubes. Physical Review Letters 2001;87:215502- 
1-4. 

[7] Yang YM, Maa JR. Boiling of suspension of solid particles in water. International 
Journal of Heat and Mass Transfer 1984;27:145-7. 

[8] Choi SUS. Enhancing thermal conductivity of fluids with nanoparticles. ASME 
FED 1995;231:99-105. 

[9] Gass V, Van der Schoot BH, de Rooij NF. Nanofluid handling by micro-flow- 
sensor based on drag force measurements. In: Proceedings of IEEE conference 
on micro electro mechanical systems. 1993. p. 167-72. 

[10] Grimm A. Powdered aluminum-containing heat transfer fluids. German patent 

DE 4131516 A1; 1993. 

[11] Lee S, Choi SUS, Li S, Eastman JA. Measuring thermal conductivity of fluids 

containing oxide nanoparticles. Journal of Heat Transfer 1999;121:280-9. 

[12] Murshed SMS, Leong KC, Yang C. Enhanced thermal conductivity of TiO2-water 
based nanofluids. International Journal of Thermal Sciences 2005;44: 
367-73. 

13] Murshed SMS, Leong KC, Yang C. Investigations of thermal conductivity and 
viscosity of nanofluids. International Journal of Thermal Sciences 2008;47: 
560-8. 

14] Yu W, France DM, Routbort JL, Choi SUS. Review and comparison of nanofluid 
thermal conductivity and heat transfer enhancements. Heat Transfer Engineer- 
ing 2008;29:432-60. 

15] Murshed SMS, Leong KC, Yang C. Thermophysical and electrokinetic 
properties of nanofluids—a critical review. Applied Thermal Engineering 
2008;28:2109-25. 

16] Murshed SMS, Leong KC, Yang C. Determination of the effective thermal dif- 
fusivity of nanofluids by the double hot-wire technique. Journal of Physics D: 
Applied Physics 2006;39:5316-22. 

17] Buongiorno, et al. A benchmark study on the thermal conductivity of nanoflu- 
ids. Journal of Applied Physics 2009;106:094312. 

18] Tong LS, Tang YS. Boiling heat transfer and two-phase flow. Taylor & Francis; 
1997. 

19] Pioro IL, Rohsenow W, Doerffer SS. Nucleate pool-boiling heat transfer. I. 
Review of parametric effects of boiling surface. International Journal of Heat 
and Mass Transfer 2004;47:5033-44. 

20] Pioro IL, Rohsenow W, Doerffer SS. Nucleate pool-boiling heat transfer. II. 
Assessment of prediction methods. International Journal of Heat and Mass 
Transfer 2004;47:5045-57. 

21] Kang MG. Experimental investigation of tube length effect on nucleate pool 
boiling heat transfer. Annals of Nuclear Energy 1998;25:295-304. 

22] Kang MG. Effect of surface roughness on pool boiling heat transfer. International 
Journal of Heat and Mass Transfer 2000;43:4073-85. 

23] Kang MG. Effects of tube inclination on pool boiling heat transfer. Nuclear 
Engineering and Design 2003;220:67-81. 

24] You SM, Kim JH, Kim KM. Effect of nanoparticles on critical heat flux of water 
in pool boiling of heat transfer. Applied Physics Letters 2003;83:3374-6. 

25] Milanova D, Kumar R. Heat transfer behavior of silica nanoparticles in pool 
boiling experiment. Journal of Heat Transfer 2008;130:042401-1-6. 

26] Das SK, Putra N, Roetzel W. Pool boiling characterization of nano-fluids. Inter- 
national Journal of Heat and Mass Transfer 2003;46:851-62. 

27] Wen D, Ding Y. Experimental investigation into the pool boiling heat trans- 
fer of aqueous based alumina nanofluids. Journal of Nanoparticle Research 
2005;7:265-74. 

28] Witharana S. Boiling of refrigerants on enhanced surfaces and boiling of 
nanofluids. PhD thesis. Royal Institute of Technology; 2003. 

29] Prakash NG, Anoop KB, Das SK. Mechanism of enhancement/deterioration of 
boiling heat transfer using stable nanoparticles suspensions over vertical tubes. 
Journal of Applied Physics 2007;102:074317-1-7. 

30] Prakash NG, Anoop KB, Sateesh G, Das SK. Effect of surface orientation on 
pool boiling heat transfer of nanoparticle suspensions. International Journal 
of Multiphase Flow 2008;34:145-60. 


[31] Soltani S, Etemad SG, Thibault J. Pool boiling heat transfer performance of 
Newtonian nanofluids. Heat and Mass Transfer 2009;45:1555-60. 

[32] Soltani S, Etemad SG, Thibault J. Pool boiling heat transfer performance of New- 
tonian nanofluids. International Communications of Heat and Mass Transfer 
2010;37:29-33. 

[33] Truong B, Hu LW, Buongiorno J, McKrell T. Modification of sandblasted plate 
heaters using nanofluids to enhance pool boiling critical heat flux. International 
Journal of Heat and Mass Transfer 2010;53:85-94. 

[34] Kim SJ, McKrell T, Buongiorno J, Hu LW. Enhancement of flow boiling crit- 
ical heat flux (CHF) in alumina/water nanofluids. Advanced Science Letters 
2009;2:100-2. 

[35] Henderson K, Park YG, Liu L, Jacobi AM. Flow-boiling heat transfer of R-134a- 
based nanofluids in a horizontal tube. International Journal of Heat and Mass 
Transfer 2010;53:944-51. 

[36] Vassallo P, Kumar R, D’Amico S. Pool boiling heat transfer experiments in 
silica-water nano-fluids. International Journal of Heat and Mass Transfer 
2004;47:407-11. 

[37] Bang IC, Chang SH. Boiling heat transfer performance and phenomena of 
Al203-water nano-fluids from a plain surface in a pool. International Journal 
of Heat and Mass Transfer 2005;48:2407-19. 

[38] Kim H, Kim J, Kim M. Experimental study on CHF characteristics of water-TiO2 
nano-fluids. Nuclear Engineering Technology 2006;39:61-8. 

[39] Kim SL, Bang IC, Buongiorno J, Hu LW. Surface wettability change during pool 
boiling of nanofluids and its effect on critical heat flux. International Journal of 
Heat and Mass Transfer 2007;50:4105-16. 

[40] Jackson J. Investigation into the pool-boiling characteristics of gold nanofluids. 
MS thesis. University of Missouri-Columbia; 2007. 

[41] Chopkar M, Das AK, Manna I, Das PK. Pool boiling heat transfer characteristics 
of ZrO2-water nanofluids from a flat surface in a pool. Heat and Mass Transfer 
2008;44:999-1004. 

[42] Liu ZH, Liao L. Sorption and agglutination phenomenon of nanofluids on a plain 
heating surface during pool boiling. International Journal of Heat and Mass 
Transfer 2008;51:2593-602. 

[43] Lv LC, Liu ZH. Boiling characteristics in small vertical tubes with closed bot- 
tom for nanofluids and nanoparticles-suspensions. Heat and Mass Transfer 
2008;45:1-9. 

[44] Kathiravan R, Kumar R, Gupta A, Chandra R, Jain PK. Pool boiling characteristics 
of carbon nanotube based nanofluids over a horizontal tube. Journal of Thermal 
Science and Engineering Applications 2009;1:022001-1-7. 

[45] Rohsenow WM. A method of correlating heat transfer data for surface boiling 
of liquids. Transactions of the ASME 1952;74:969-76. 

[46] Zuber N. Hydrodynamic aspects of boiling heat transfer. Physics and Mathe- 
matics, AEC Report No. AECU-4439; 1959. 

[47] Pioro IL. Experimental evaluation of constants for the Rohsenow pool boil- 
ing correlation. International Journal of Heat and Mass Transfer 1999;42: 
2003-13. 

[48] Murshed SMS, Milanova D, Kumar R. An experimental study of surface tension- 
dependent pool boiling characteristics of carbon nanotubes-nanofluids. 
In: Proceedings of 7th international ASME conference on nanochannels, 
microchannels and minichannels. 2009. 

[49] Milanova D, Kumar R. Functionalized single walled and double walled carbon 
nanotubes for thermal enhancement. In: Proceedings of the ASME international 
mechanical engineering congress and exposition. 2007. 

[50] Kim H, Kim J, Kim MH. Effect of nanoparticles on CHF enhancement in 
pool boiling of nano-fluids. International Journal of Heat and Mass Transfer 
2006;49:5070-4. 

[51] Wasan DT, Nikolov AD. Spreading of nanofluids on solids. Nature 
2003;423:156-9. 

[52] Duursma G, Sefiane K, Kennedy A. Experimental studies of nanofluid droplets 
in spray cooling. Heat Transfer Engineering 2009;30:1108-20. 

[53] Shen J, Liburdy JA, Pence DV, Narayanan V. Droplet impingement dynamics: 
effect of surface temperature during boiling and non-boiling conditions. Journal 
of Physics: Condensed Matter 2009;21:464133-1-14. 

[54] Murshed SMS, Nieto de Castro CA. Spreading characteristics of nanofluid 
droplets impacting onto a solid surface. Journal of Nanoscience and Nanotech- 
nology 2011;11:1-7. 

[55] Pak BC, Cho YI. Hydrodynamic and heat transfer study of dispersed flu- 
ids with submicron metallic oxide particles. Experimental Heat Transfer 
1998;11:151-70. 

[56] Eastman JA, Choi SUS, Li S, Soyez G, Thompson LJ, Dimelfi RJ. Novel 
thermal properties of nanostructured materials. Materials Science Forum 
1999;312-314:629-34. 

[57] Xuan Y, Li Q. Investigation on convective heat transfer and flow features of 
nanofluids. Journal of Heat Transfer 2003;125:151-5. 

[58] Wen D, Ding Y. Experimental investigation into convective heat transfer of 
nanofluids at the entrance region under laminar flow conditions. International 
Journal of Heat and Mass Transfer 2004;47:5181-8. 

[59] Heris SZ, Etemad SG, Esfahany MS. Experimental investigation of oxide nanoflu- 
ids under laminar flow convective heat transfer. International Communications 
of Heat and Mass Transfer 2006;33:529-35. 

[60] Duangthongsuk W, Wongwises S. Heat transfer enhancement and pressure 
drop characteristics of TiO2-water nanofluid in a double-tube counter flow heat 
exchanger. International Journal of Heat and Mass Transfer 2009;52:2059-67. 

[61] Faulkner D, Rector DR, Davison JJ, Shekarriz R. Enhanced heat transfer through 
the use of nanofluids in forced convection. In: Proceedings of ASME heat trans- 
fer div. 2004. p. 219-24. 


2354 S.M.S. Murshed et al. / Renewable and Sustainable Energy Reviews 15 (2011) 2342-2354 


[62] Jung JY, Oh HS, Kwak HY. Forced convective heat transfer of nanofluids in 
microchannels. In: Proceedings of the ASME international mechanical engi- 
neering congress and exposition. 2006. 

[63] Ding Y, Alias H, Wen D, Williams AR. Heat transfer of aqueous suspen- 
sions of carbon nanotubes. International Journal of Heat and Mass Transfer 
2006;49:240-50. 

[64] Yang Y, Zhang ZG, Grulke EA, Anderson WB, Wu G. Heat transfer properties 
of nanoparticle-in-fluid dispersions (nanofluids) in laminar flow. International 
Journal of Heat and Mass Transfer 2005;48:1107-16. 

[65] Lai WY, Duculescu B, Phelan PE, Prasher RS. Convective heat transfer with 
nanofluids in a single 1.02-mm tube. In: Proceedings of the ASME international 
mechanical engineering congress and exposition. 2006. 

[66] Williams W, Buongiorno J, Hu LW. Experimental investigation of turbulent con- 
vective heat transfer and pressure loss of alumina/water and zirconia/water 
nanoparticle colloids (nanofluids) in horizontal tubes. Journal of Heat Transfer 
2008;130:1-7. 

[67] Hwang KS, Jang SP, Choi SUS. Flow and convective heat transfer characteris- 
tics of water-based Al203 nanofluids in fully developed laminar flow regime. 
International Journal of Heat and Mass Transfer 2009;52:193-9. 

[68] Xie H, Li Y, Yu W. Intriguingly high convective heat transfer enhancement of 
nanofluid coolants in laminar flows. Physics Letters A 2010;374:2566-8. 

[69] Amrollahi A, Rashidi AM, Lotfi R, Meibodi ME, Kashefi K. Convection heat trans- 
fer of functionalized MWNT in aqueous fluids in laminar and turbulent flow at 
the entrance region. International Communications in Heat and Mass Transfer 
2010;37:717-23. 

[70] Kulkarni DP, Namburu PK, Bargar HE, Das DK. Convective heat transfer and fluid 
dynamic characteristics of SiO2-ethylene glycol/water nanofluid. Heat Transfer 
Engineering 2008;29:1027-35. 

[71] Daungthongsuk W, Wongwises S. A critical review of convective heat transfer 
of nanofluids. Renewable and Sustainable Energy Reviews 2007;11:797-817. 

[72] Kakag S, Pramuanjaroenkij A. Review of convective heat transfer enhance- 
ment with nanofluids. International Journal of Heat and Mass Transfer 
2009;52:3187-96. 

[73] Putra N, Roetzel W, Das SK. Natural convection of nanofluids. Heat and Mass 
Transfer 2003;39:775-84. 

[74] Khanafer K, Vafai K, Lightstone M. Buoyancy-driven heat transfer enhancement 
in a two-dimensional enclosure utilizing nanofluids. International Journal of 
Heat and Mass Transfer 2003;46:3639-53. 

[75] Wen D, Ding Y. Formulation of nanofluids for natural convective heat transfer 
applications. International Journal of Heat and Fluid Flow 2005;26:855-64. 

[76] Buongiorno J. Convective transport in nanofluids. Journal of Heat Transfer 
2006; 128:240-50. 

[77] Bejan A. Convection heat transfer. John Wily & Sons, Inc.; 2004. 


[78] Dittus FW, Boelter LMK. Heat transfer in automobile radiators of the tubular 
type. University of California Publications on Engineering 1930;2:443-61. 

[79] Maiga EEB, Palm SJ, Nguyen CT, Roy G, Galanis N. Heat transfer enhancement 
by using nanofluids in forced convection flows. International Journal of Heat 
and Fluid Flow 2005;26:530-46. 

[80] Roy G, Nguyen CT, Lajoie PR. Numerical investigation of laminar flow and heat 
transfer in a radial flow cooling system with the use of nanofluids. Superlattices 
and Microstructures 2004;35:497-511. 

[81] Xuan Y, Yu K, Li Q. Investigation on flow and heat transfer of nanofluids by the 
thermal Lattice Boltzmann model. Progress in Computational Fluid Dynamics 
2005;5:13-9. 

[82] Palm SJ, Roy G, Nguyen CT. Heat transfer enhancement with the use of 
nanofluids in radial flow cooling systems considering temperature-dependent 
properties. Applied Thermal Engineering 2006;26:2209-18. 

[83] Jou RY, Tzeng SC. Numerical research of nature convective heat transfer 
enhancement filled with nanofluids in rectangular enclosures. International 
Communications in Heat and Mass Transfer 2006;33:727-36. 

[84] Jang SP, Choi SUS. Cooling performance of a microchannel heat sink with 
nanofluids. Applied Thermal Engineering 2006;26:2457-63. 

[85] Heris SZ, Esfahany MN, Etemad G. Numerical investigation of nanofluid laminar 
convective heat transfer through a circular tube. Numerical Heat Transfer Part 
A: Applications 2007;52:1043-58. 

[86] MirmasoumiS, Behzadmehr A. Effect of nanoparticles mean diameter on mixed 
convection heat transfer of a nanofluid in a horizontal tube. International Jour- 
nal of Heat and Fluid Flow 2008;29:557-66. 

[87] Ho CJ, Chen MW, Li ZW. Numerical simulation of natural convection of 
nanofluid in a square enclosure: effects due to uncertainties of viscosity 
and thermal conductivity. International Journal of Heat and Mass Transfer 
2008;51:4506-16. 

[88] Oztop HF, Abu-Nada E. Numerical study of natural convection in partially 
heated rectangular enclosures filled with nanofluids. International Journal of 
Heat and Fluid Flow 2008;29:1326-36. 

[89] Vajjha RS, Das DK, Namburu PK. Numerical study of fluid dynamic and heat 
transfer performance of Al203 and CuO nanofluids in the flat tubes of a radiator. 
International Journal of Heat and Fluid Flow 2010;31:613-21. 

[90] Murshed SMS, Leong KC, Yang C, Nguyen NT. Convective heat transfer 
characteristics of aqueous TiO nanofluids under laminar flow conditions. 
International Journal of Nanoscience 2008;7:325-31. 

[91] Hamilton RL, Crosser OK. Thermal conductivity of heterogeneous two 
component systems. Industrial and Engineering Chemistry Fundamentals 
1962;1:187-91. 

[92] Batchelor GK. The effect of Brownian motion on the bulk stress in a suspension 
of spherical particles. Journal of Fluid Mechanics 1977;83:97-117. 


